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Abstract. We present some recent theory and simulation results addressing 
the phenomena of colloidal gelation at both high and low volume fractions, 
in the presence of short-range attractive interactions. We discuss the ability 
of mode-coupling theory and its adaptations to address situations with strong 
heterogeneity in density and/or dynamics. We include a discussion of the effect 
of attractions on the shear-thinning and yield behaviour under flow. 
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1. Introduction 

Recent studies of mode-coupling theory (MCT) have predicted various kinetic arrest 
scenarios for colloids with short-range attractions [1-3]. The behaviour at the arrest 
transition has been analysed in considerable detail as a function of the volume fraction 
4> of the colloids, their attraction range 5 (in units of the particle radius a) and the 
attraction e (in units of fcgT) [4,5]. This approach, which builds on a successful 
theory of the glass transition for hard spheres, suggests that the underlying 'ideal glass 
transition', found within MCT, could provide a universal mechanism for homogeneous 
gelation in which the arrested state is viewed as an attraction-driven glass [6, 7]. 

Here we do not attempt to mediate between those who support and those who 
oppose the basic philosophy of MCT (see [8, 9] for discussions). Instead, in Sections 2- 
4, we outline three areas of recent work that examine the ability of MCT to deal with 
heterogeneous systems. In Section 2 we consider attractive colloids at low density, 
and consider kinetic arrest by routes that combine aspects of irreversible aggregation, 
gelation, and phase separation. The phenomenology that we describe takes MCT as 
a reliable theory of arrest at high densities only, and attempts to extend the same 
picture to the much lower colloid volume fractions at which gelation actually occurs in 
systems with strong short-range attractions. The resulting gels are very heterogeneous 
and our approach is to preserve what we can of the MCT theory while allowing for 
heterogeneity in density at scales either shorter or longer than those at which we apply 
the MCT. In the first case (shorter scale heterogeneity) one has a picture involving the 
MCT-like arrest of clusters; in the second (longer scale) one has a phase-separation 
driven morphology in which one phase arrests into a dense gel or attractive glass. 

In Section 3 we turn to gelation at higher particle densities, where unmodified 
MCT can reasonably hope to succeed, at least in the absence of phase separation. 
Here we carefully analyse simulation data for attractive colloids, suppressing phase 
separation by including a weak, long-range repulsion. We find that MCT gives 
an excellent account of averaged dynamical quantities but that this disguises an 
underlying dynamics which is much richer, and which in fact shows many hallmarks of 
dynamical heterogeneity (DH): populations of fast and slow-moving particles coexist. 
This heterogeneity appears to be closely related to static heterogeneity of density 
which is present for attractive colloids even at quite high volume fractions. This form 
of DH may thus be different from that found in hard spheres [10] where the density 
fluctuations are smaller. The new type of DH may be enhanced by the long-range 
repulsion we have added, but even if that turns out to be true, the results are of 
strong interest. They show firstly that when experimental data agrees with MCT 
this does not exclude the possibility of DH; likewise that observation of DH does not 
exclude the possibility that MCT remains predictive for averaged quantities such as 
the dynamic structure factor on the fluid side of the transition. 

In Section 4 we describe some related recent work on colloids under shear, and give 
a preliminary account of calculations that address the effect of short range attractions 
on the rheological behaviour. This offers nontrivial predictions, in particular for the 
variation of yield stress and shear thinning behaviour with the range of attraction S. 
In Section 5 we offer some concluding remarks. 
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2. Gelation at low density 

2.1. MCT, phase separation, and aggregation 

For attractive colloids, MCT appears to be very useful at predicting averaged 
dynamical quantities like the dynamic structure factor S(q,t), so long as the volume 
fraction <j> of colloid is high enough [7, 11, 12]. In the following we refer to the resulting 
spatially homogenous gel as type I. In practice, for moderately short relative range 
of attraction 5 > 10~ 2 a metastable gas-liquid phase separation [13] interferes with 
this kinetic arrest scenario. On quenching to create a gel, phase separation is likely 
to intervene for all colloid volume fractions below the intersection of the binodal with 
the MCT arrest line. 

The 'hidden' binodal, when present, dominates over the much slower 
crystallization route to the thermodynamic ground state, so that it is liquid-gas 
separation not crystallization that can interfere more strongly with glass formation. 
Moreover, due to the kinetic (as opposed to thermodynamic) mechanism underlying 
the arrest transition, this potentially complex interplay sensitively depends on the 
quench rate relative to the typical time scale for phase separation. To a first 
approximation, the latter is given by the Smoluchowski time, which is the time scale 
for doublets to form by colloidal collisions. Depending on whether the quench is fast 
or slow with respect to this natural time scale, we expect transient gels to develop 
along one of two different routes, in the following referred to as type II and type III, 
respectively. Type-II gels are homogeneous on short scales but strongly heterogeneous 
at the mesoscale: they result when the characteristic coarsening textures produced 
by phase separation get 'frozen in' during the coarsening process, as a result of 
an MCT-like arrest of one of the two phases [14]. Type-Ill gels are, in contrast, 
assemblies of long-lived nonequilibrium structures locally resembling those obtained 
from irreversible cluster aggregation [15] and thus heterogeneous also on short scales. 
To distinguish from fully irreversible aggregates, in which the bonds formed are 
permanent, the type III process is sometimes called 'weak gelation'. 

During the process leading to type-II gels, the system remains in (or sufficiently 
close to) local equilibrium, so that gelation by this route represents a relatively 
straightforward combination of phase separation and MCT. However, additional 
concepts are needed for the type-Ill scenario. Some of us have recently proposed 
a schematic description (called cluster-MCT or CMCT) [16] for weak gelation in such 
suspensions, in terms of an effective theory for a fluid of 'renormalised particles'. These 
represent coarse-grained clusters and their interactions are analysed, in an MCT-like 
fashion, to predict the onset of global kinetic arrest. The underlying view of gelation 
as a double ergodicity breaking (on the monomer scale and on the cluster scale) seems 
to be supported by recent numerical simulations [17]. 

2.2. Schematic nonequilibrium phase diagram 

A schematic phase diagram of the different nonequilibrium behaviours predicted for 
hard sphere colloids with short-range attractions is provided in Fig. 1. 

Slow quench scenario: We concentrate first on this conceptually simple case, 
which results in type-II gels. Consider a slow quench along the quench path indicated 
by the arrow at volume fraction <j) a . Upon crossing the (metastable) liquid-gas binodal, 
the fluid starts to decompose into a minority phase consisting of colloidal 'liquid' 
domains (rich in particles), and a majority 'gas' phase consisting of domains where 
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Figure 1. Schematic [cj> — e) cut through the phase diagram of adhesive hard 
spheres. The interaction of gas-liquid phase separation with MCT kinetic arrest 
gives rise to a variety of metastable gel phases. At high colloid volume fractions <fi, 
where MCT can be applied directly, homogeneous gels form (region I). At lower 
<f> phase separation may create macroscopically heterogeneous gels (region II) and 
(homogeneous) gel beads if the quench is slow, while more tenuous gels (III) and 
(ramified) clusters consisting of non— equilibrium particle aggregates result from 
rapid quenches. 

colloids are scarce. These domains undergo slow coarsening, while their compositions 
continue to evolve along the binodal in response to the ongoing quench. During this 
process, the two phases both remain in local coexistence on the binodal and move apart 
until the compositions (pi and tf2 are reached. Note that, although the quench is slow, 
we assume that macroscopic phase separation is even slower, so that the morphology 
is one of mesoscopic domains. Once the colloid volume fraction within the minority 
phase reaches ip2 (see Fig. 1), which is where the binodal cuts the MCT line, the denser 
domains undergo kinetic arrest: the characteristic time scales for further structural 
evolution suddenly increase dramatically. From the lever rule, the volume fraction of 
space occupied by the resulting amorphous solid phase is = ((f) — tp 1 )/(ip 2 — (fi). 

Depending on the mesoscale domain structure that has accompanied the 
preceding phase separation, this volume fraction $ may be more or less than the value 
$ c required for the minority phase to percolate. Only if it does percolate do we have 
a macroscopic solid or gel: note that gelation, in the sense of a finite elastic modulus, 
requires percolation of an arrested phase, or equivalently arrest of a percolated one. 
In particular, existence in the system of a percolating network of bonds, each of which 
is transient, is not sufficient to create a modulus. 
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In terms of the colloid volume fraction <fi, we obtain a percolation threshold (j> c 
defined by <fr((j) c ) — 3> c . For a volume fraction fa > 4> c as depicted in the figure, 
a connected gel, of finite modulus, is the predicted result of a quench. By varying 
the colloid volume fraction and the quench kinetics, different gel morphologies should 
be realizable in a relatively well-controlled manner. Indeed, corresponding recipes 
are routinely applied in the industrial processing of colloidal and polymeric gels for 
optimizing rheological properties (see e.g. [18,19]). Upon decreasing <fi, the gel gets 
more tencous and eventually, for <f> < <f> c (e.g., <f> — <f> a in Fig. 1) this disintegrates into 
a fluid of slightly sticky gel beads. (At still lower densities, there is a slight shift in 
the binodal curve because the coexisting dense phase cannot be denser than (j> 2 .) 

For moderately deep quenches, the mutual attraction between the colloids, though 
strong enough to solidify dense colloidal drops (forming the aforementioned beads), 
is insufficient to permanently bind beads that come into contact. This is because 
the beads are now made of colloidal glass, and cannot adapt their shapes to allow 
coalescence on the Smoluchowski timescale, even though this would reduce their 
surface energy significantly. Two such beads of radius R, treated as effective particles, 
have a far smaller relative range of attraction than the primary particles (S c s ~ S/R), 
which for roughened surfaces may be only partly outweighed by its increased depth 
e c ff t- For deeper quenches, the residual interaction between beads may become strong 
enough to cause a second round of phase separation and/or gelation at the bead level; 
the whole argument can be iterated to describe that case. 

Ignoring this last effect for the moment, the gel region in the phase diagram 
is bounded by a section of liquid-gas tieline emanating from the intersection point 
(at ip 2 ) of the binodal with the MCT arrest line, and a line parallel to the quench 
route, emanating from this tieline at the volume fraction <p — <Pc = fi + $ C (V2 — 
ipi) corresponding to percolation. The transition across the first boundary is 
'temperature-driven' (or more accurately ^-driven) and governed by the MCT arrest 
of the dense phase; quenching by this route, some features of the MCT transition 
(which applies directly only for (f> > ip?) are expected to be detectable dynamically at 
much lower concentrations [16]. In contrast, the 'pressure-driven' (or more accurately, 
0-driven) transition across the remaining section of the gel boundary should fall into 
the percolation universality class [20]. The possible secondary gelation of arrested 
beads at deep quenches complicates this picture somewhat; we have neglected it in 
the figure and do not pursue this here. 

Rapid quench scenario: Rapid quenches are ones in which the system finds itself 
under conditions of strong, quasi-irreversible aggregation, so that bonds created in 
collisions between particles are very long lived. In this case one expects at low volume 
fractions a nontrivial episode of structure formation, akin to irreversible diffusion- 
limited aggregation, on a timescale fast compared to other processes. Morphologically, 
rapid quenches should differ from the slow quenches described above, with a more 
ramified local structure and a correspondingly larger region in which an arrested phase 
percolates; the gel phase thus extends to lower volume fractions than in a slow quench. 

The resulting type-Ill gels are conceptually distinct from the rest, and may have 
some intriguing properties. Their analysis is more difficult, as can be appreciated by 
considering a deep rapid quench at low volume fractions. This is the case considered 

J Note that for ideal spheres £ e g would typically increase linearly with size and thus over-compensate 
the effect of the reduction in the effective range S c g . However, the bead surfaces are expected to be 
roughened due to particle deposition in the ongoing quench, so that estimating the effective attraction 
is a subtle task left for future work. 
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in the CMCT theory of Ref.[16]. The basic idea of CMCT is to allow for strong, 
quasi-irreversible bonding at short length scales by applying MCT, not to primary 
particles, but to fractal aggregates created by such bonding. This is quite similar to 
the idea already introduced above, of iterating MCT to describe the possible gelation 
and arrest of dense beads. However, the calculation of the effective particle parameters 
is even more subtle and, in particular, allows for an entropic decrease in the effective 
attraction between clusters. This reflects the multiplicity of internal bonds at which 
two clusters, having joined to form a larger one, can now be broken apart to recover 
two clusters of similar size and shape to the original pair [16]. Another potentially 
important feature is the ability of CMCT to allow for a buildup of long-range repulsions 
due, for example, to a very slight Coulomb repulsion that may be present in many 
colloids, even in organic solvents [21]. This buildup also lowers the effective attraction 
strength as the clusters get larger. 

The effect of the scale-dependent effective bond strength, combined with a scale- 
dependent range, means that in CMCT there is a tendency, as aggregation proceeds, to 
move away from the attraction-driven arrest scenario and towards a more conventional 
repulsion-driven glass. Setting aside phase separation effects (for now), the key issue 
in gelation is then whether the effective attractions become small while the volume 
fraction of clusters is still fairly low, or whether, by the time they stop aggregating, the 
clusters are dense enough to be arrested anyway by repulsive caging. In the first case 
one predicts a semi-ergodic phase comprising a fluid of clusters [22]; in the second, an 
arrested cluster phase which is nonergodic at all scales and thus a gel [16]. 

For type-Ill gels the interplay with phase separation is quite complicated; roughly 
speaking it follows the lines already developed above for slow quenches, but with 
CMCT-like arrest replacing the standard MCT arrest throughout the discussion. 
However the phase separation itself is also perturbed by what is happening at the 
cluster scale, which is a further complication. (In mitigation, for many systems the 
dense phase will be dense enough for the CMCT and MCT predictions to nearly 
coincide through much of the phase diagram anyway.) Within the CMCT picture 
there can be two distinct forms of cluster phase: one in which there is no tendency to 
macroscopic phase separation but clusters stop growing due to the buildup of repulsion 
(Coulombic or entropic); another in which phases have separated but the dense 
arrested phase does not percolate [16]. The latter includes the gel beads described 
previously as a limiting case [14]. 

CMCT is based on a simplifying assumption that the timescale for internal 
reconstruction of clusters is slow compared to the timescale for realizing a state of 
repulsion-driven arrest at larger scales. This may be safe for Coulombic stabilization, 
but needs careful further test, against both experiment and simulation, in the case 
where the bond strength is effectively reduced by the entropy associated with bond 
breaking internal to a cluster [23] . 

2.3. Combining MCT with phase separation: further remarks 

By translating knowledge of the interaction potential into a combined topographic map 
of the MCT transition surface and the gas-liquid coexistence, we were able to present 
above and in [16] (see also [14]) a guide for disentangling some of the experimentally 
observed complex phenomenology of attracting colloids. 

We now take a somewhat broader perspective in which, on top of the hard sphere 
repulsions, the pair potential can be represented as the sum of a strong short-range 
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Figure 2. Static structure factor S q for the AHSS [24] with volume fraction 
<p = 0.3 and adhesiveness parameter r = 0.08 and r = 0.2, respectively 
(the former being closer to the spinodal and showing the more pronounced 
oscillations), calculated by a random phase approximation. Three cases of 
competing interactions are distinguished: The bare AHSS (dashed lines), the 
AHSS with an additional repulsive Yukawa barrier of height 0.1 and decay length 
2.5 (solid lines), and with an attractive tail of the same range but with an 
amplitude —0.004 (dotted line, only for r = 0.2 since r = 0.08 lies within the 
spinodal region). 

attraction and an additional weak long-range interaction that may either be attractive 
of repulsive. Numerous examples — from microemulsions, through block copolymers, 
to supercooled water — show that such competing interactions can give rise to complex 
phase behavior. Even the adhesive hard-sphere system (AHSS) of Baxter [24] is now 
known to have two competing crystalline phases, two glasses, and a nontrivial re-entry 
between them. This stems from a competition between the hard-core repulsion and 
the short range attraction [6, 25]. Adding further long-range interactions can certainly 
complicate this further. 

Insofar as it destabilizes the equilibrium liquid, the origin of this complexity 
is manifest in the static structure factor S q [26]; this is also the input for MCT 
calculations. In Fig. 2, S q is shown for an AHSS with additional long range interaction. 
The low— q tail of the structure factor shows the characteristic upturn found near any 
spinodal. A sharper upturn is seen (well before the spinodal is reached) if a very 
weak but longer ranged attraction is added (dotted). In contrast, adding a repulsive 
barrier can energetically prohibit macroscopic phase separation as evidenced by the 
suppression of the spinodal divergence and the shifting of the new relative maximum 
in S q to finite wave vectors q = q c (solid lines). Near the former spinodal, the new 
peak grows in height without bound, thus heralding microphase separation, in which 
the homogeneous fluid decomposes into finite, self-limiting domains of liquid-like and 
gas-like character. The domain size (~ q^ 1 ) is controlled by (and roughly scales as) 
the range of the extra repulsion. 

To predict the non-equilibrium behaviour in the presence of the competing 
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interactions, note first that thermodynamics and kinetic arrest are sensitive to 
different features of the pair potential. Gas-liquid phase separation can be induced 
by a weak long-ranged attraction, as familiar from van der Waals' theory; and 
similarly, even finite repulsions may cause (Wigner-)crystallization if sufficiently long- 
ranged. Attraction-driven arrest, on the other hand, relies on kinetic trapping of a 
thermodynamically stable fluid by short-ranged forces (S <C 1), and according to MCT 
the arrest is triggered by the high Fourier modes qa ~ S^ 1 » 1, as encoded in the 
direct correlation function c q oc 1 — Sg 1 [1]. This contrasts with the crowding of the 
nearest-neighbour shell visible in the first peak of S q near qa ~ n, which leads to the 
"cage effect" for repulsive hard sphere glasses. 

Recently there has been much debate whether the long range repulsion itself 
[27], or the frustration resulting from the competing interactions [28], may be a 
new driving mechanism for the formation of more exotic arrested states. Indeed, 
apart from the repulsive (and attractive) hard-core glasses familiar in colloids, MCT 
allows for another transition scenario triggered by long range repulsions that give 
rise to a dilute arrested state, which was interpreted as a Wigner glass [29]. Even 
liquid-crystalline domains of asymmetric particles have been predicted to undergo 
kinetic arrest [30]. Along somewhat different routes similar conclusions were reached 
about so-called stripe glasses in microemulsions and related systems [28]. Several of 
the mentioned mechanisms (attractive/repulsive 'hard-core' interactions of arrested 
spinodal textures; a glass transition triggered by a microphase peak in S q ; the 
Wigner glass) may be relevant to the possible gel phases of Coulomb-stabilised cluster 
fluids [21,22]. More generally we therefore expect that the interplay of MCT with 
microphase separation could lead to a phenomenology at least as rich and interesting 
as the one elaborated above for the case of bulk phase separation. 

3. Simulation of dense attractive colloids 

We now turn to our second theme, which is to test the predictions of MCT for 
dense attractive colloids, in a regime of concentration and interactions where these 
predictions could be reasonably expected to work. Thus we can bypass the various 
complications connected with phase separation and aggregation that arise at low 
density, as were considered in Section 2. 

3.1. Simulation details 

Newtonian dynamics simulations were performed to test the theoretical predictions 
on the gel transition in a system whose pair potential mimics the depletion attraction 
found in colloid polymer mixtures [6]. One thousand polydisperse particles were 
considered in the canonical ensemble. The core-core repulsion between particles is 
given by V sc (r) — ksT (r /a^)^ 36 , where ai2 = ai + a2, with a\ and <22 the radii of 
the particles. A flat distribution of radii was used to prevent crystallisation, with a 
(half-)width equal to one tenth of mean radius, A = 0.1a. The interaction between the 
colloidal spheres is given by the Asakura-Oosawa interaction potential, which considers 
the depletion of ideal polymers [31], corrected to take into account the polydispersity 
of the colloids [32]. The total potential (T^ot = Vao + V sc ) was corrected close to a\ 2 
to ensure that the minimum of the total potential is at a\i [33]. 

A long-range repulsive barrier was added to the interaction potential in order 
to prevent liquid-gas separation at high attraction strength. When this barrier is 
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Figure 3. Long time limit of the demixing order parameter, i/>4 as a function of 
the polymer volume fraction, <j> p . Inset: Evolution of ip4 with time for different 
<b v as labelled. 



allowed for, the potential is attractive at short distances, and slightly repulsive at 
longer distances. The range of the attraction is given by the size of the polymers, £, 
and the strength at contact, r = an, is proportional to the polymer volume fraction 
(pp. The height of the repulsive barrier is lksT whereas arrest occurs for attraction 
strengths above ~ 7ksT [33]. While the remarks of Section 2.3 serve as a warning 
that this kind of barrier could influence the results, our aim here is merely to avoid 
the liquid-gas separation that would otherwise complicate matters considerably. The 
barrier is thus chosen as small as is consistent with achieving this. 

In our simulations, lengths are measured in units of the mean particle radius, a, 
and time in units of ^/4a 2 /3w 2 , with v the thermal velocity, which is set to \/4/3. 
The equations of motion were integrated using the velocity- Verlet algorithm, with a 
time step of 0.0025. The volume fraction of colloids is fixed at <f> = 0.40; £ = 0.1a; and 
the attraction strength is parameterised by the polymer volume fraction, <p p . 

3.2. Suppression of phase separation 

Separation of the system into two phases of different density was monitored by a 
demixing order parameter, -04- The system is divided in 4 3 boxes of equal size, and 
the deviation of the density in every box with respect to the average is measured by 
this parameter; an homogeneous system shows ip4 close to zero and a phase-separated 
system has a much bigger value. With this parameter, the isochore (f> = 0.40 was first 
studied for the system without a long range barrier, in order to see any effect gelation 
might have on liquid-gas separation (and, potentially, vice versa). 

The inset to figure 3 shows the temporal evolution of V4 for different polymer 
fractions. At low <f> p , the system is homogeneous, i.e., the attraction strength is too 
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low to induce liquid-gas separation, and for <j> p = 0.35, the system demixes, as shown 
by the increase of 1^4 ■ At even higher attraction strengths, the system does not phase 
separate into denser and more dilute phases, as expected, but becomes again more 
homogeneous. The long time limit of the separation parameter, which is plotted 
in figure 3 as a function of (f> p , captures this scenario. The liquid-gas transition 
is found to take place at <f> p — 0.30, in agreement with previous simulations [34]. 
The V4 parameter does not increase monotonically, as expected from equilibrium 
thermodynamics, but reaches a maximum and decreases. A plateau at ^4 ~ 0.08 
is observed above <p p = 0.60, indicating a quasi-homogeneous system. We may thus 
conclude that some mechanism must be present that hinders liquid gas separation 
at high attraction strength and finally prevents it. This mechanism is gelation, i.e., 
attraction-driven arrest. 

3. 3. Results for averaged quantities 

To analyze the arrest mechanism, the long-range repulsive barrier was restored to the 
interaction potential, thus forbidding liquid-gas separation and allowing the gelation 
to be probed in its own right. With the barrier, the system is macroscopically 
homogeneous at all polymer fractions studied, although it presents detectable voids 
and 'tunnels' when viewed from certain angles, and also shows a low-angle peak in 
the structure factor at qa ~ 1 [33]. However this is always lower than the primary 
near-neighbour peak at qa ~ ir. 

The dynamics of the system close to gelation is studied by means of the self part 
of the density-density correlation function, i.e. 3>g(i) = (cxp{iq. [vj{t) — i-j(O)]}), 
where the brackets indicate averaging over particle j and time origin, and q is the 
wavevector. Figure 4 shows the density correlator for increasing polymer fraction for 
qa = 6.9, the second (non-microphase) peak in the structure factor. The two step 
decay in the correlation functions is similar to the behaviour of states approaching the 
glass transition in the Lennard- Jones system (LJS) or in hard spheres (HS) [35,36]. 
Furthermore, the decay from the plateau can be rescaled at long times for all the 
states presented (with deviations for the state closest to the transition, as seen also 
in other simulations [35, 36]) and the master decay can be fitted using the Kohlrausch 
form, &g(t) = A q exp(— (t/r/f)^), a signature of MCT-like nonergodicity transitions. 
Similar scalings are obtained at all the wavevectors studied (not presented here). 

Using the predictions from MCT, the early decay from the plateau is correctly 
described by the von Schwcidler law: <f> s q (t) = f s q - h < q 1) (t/T) b + h ( q \t/T) 2b + 0(t 3h ), 

where /| is the non-ergodicity parameter and h^P and are amplitudes. All three 
of these are state- independent (specifically, independent of 4> p ), whereas t is a time 
scale which carries state-dependent information, increasing as the glass transition 
is approached. The von Schweidler expression correctly describes our correlation 
functions for all states and wavevectors, and the results from the fittings show that 
f* in this attractive case is much bigger than the nonergodicity parameters found in 
HS or LJS. This fact indicates that the localization length is much shorter than in 
those cases, showing that the driving mechanism for this transition is the formation of 
long-lived bonds between particles. Moreover, the von Schweidler exponent, b, from 
the fittings also differs from the HS or LJS values, yielding b = 0.37 [33], in agreement 
with the predictions from MCT [2]. 

For the time scales r q of the a decay, defined by & q {T q ) = f q /e, MCT predicts 
a power law divergence, with an exponent 7 related to b. Figure 5 presents the time 
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Figure 4. Density-density correlation functions for increasing polymer fraction, 
rescaled at long times. From left to right <f> p = 0.425, 0.42, 0.415, 0.41, 0.40.0.39 
and 0.375, and Kohlrausch fitting to the a-decay (dashed curve). Inset: Unsealed 
correlation functions for the same states, from right to left. The hard-sphere case, 
4> P = is added (left curve). 

scales at different wavevectors on logarithmic axes. The power-law divergence is very 
clear, and the exponents agree with the value expected from 6, which is 7 = 3.1. The 
gel point, (ftp , also fitted in the analysis, is found at <j>p = 0.4265 for this volume 
fraction, <f> = 0.40. It is interesting to note that the competition between gelation 
and liquid-gas separation at this polymer fraction was already apparent in our study 
of the system without the long range barrier, although gelation fully impedes phase 
separation only well above this value of <j> p . 

MCT also predicts a power-law decay, with the same exponent 7, for the self- 
diffusion coefficient, Dq. This coefficient is determined from the long time behaviour 
of the mean squared displacement, (5r 2 ) = 6Dot. Fixing the gel point to the value 
reported above, the fitted exponent in this case is 7 = 1.23, quite different from 
the value obtained above for the divergence of the time scale. Such differences are 
obtained in the analysis of other model systems, LJS or HS, and have been attributed 
to the presence of dynamical heterogeneities in the system. However in our work the 
difference between the exponents determined from the time scale and the diffusion 
coefficient is larger than in these other cases. This suggests a possible stronger role 
for dynamic heterogeneities in the presence of short-range attractions. 

3. 4- Dynamical heterogeneity 

The dynamical heterogeneities of a state can be studied by analysing the distribution of 
the squared displacement of particles measured between some arbitrary time t = and 
a later time t — t* . A homogeneous fluid should present a single peaked distribution, 
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Figure 5. Wavevector dependent time scale, r q vs. <f>^ — <f> p for different 
wavevectors. From top to bottom qa = 3.9, 6.9, 9.9, 15, 20, 25, 30. The lines are 
power law fittings. Inset: Self diffusion coefficient, D s , vs. <f>p — <f> p and power 
law fitting. 

its width depending on t* and Dq. Figure 6 presents this distribution for different 
states approaching the gel transition, where t* has been chosen so that (5r 2 (t*)) = 10a 2 
for all states. At low polymer fractions, the system is indeed homogeneous and the 
distribution of squared displacements agrees with the theoretical expectation for a 
system of Brownian particles. However, at higher attraction strengths, the system 
becomes more and more heterogeneous and the distribution more bimodal. 

For polymer volume fractions above <f> p — 0.41 two peaks are clearly observed 
in the distribution, showing one population of particles that move much less than 
a particle radius and another population of very mobile particles. (These of course 
combine such that the average squared displacement is (8r 2 ) = 10a 2 , since this was 
used to define t* .) In fact, these two populations of particles can be distinguished at all 
times. The exchange between them is detectable, but very slow, allowing an analysis 
of the system to proceed as if it were composed of two entirely distinct populations: 
'fast' particles and 'slow' ones [37]. 

The nontrivial structure of the incipient gel (with voids and tunnels) makes 
possible the existence of a subset of particles fully integrated into a percolating gel- like 
structure, each with a large number of neighbours, while also allowing another subset 
of particles to be present at the surface of the same structure and thus to have only a 
few neighbours. The slow particles form a stiff structure which is very slow to relax, 
whereas the fast particles surround this structure, with a high mobility. Since there is 
a relatively large amount of free space in the system, the fast particles move freely and 
not in clusters or in string like motions. The fast particles attach from time to time 
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Figure 6. Distribution of squared displacement for (f> p = 0.30 (gray dotted line), 
<j> p = 0.40 (gray dashed line), (f> p = 0.41 (gray solid line), (f> p = 0.415 (black 
dotted line), <p p = 0.42 (black dashed line), and <fi p = 0.425 (black solid line). 
In all cases (<5r 2 ) = 10a 2 . The thin black line is the theoretical distribution for 
a system of non-interacting Brownian particles. Inset: Number of slow particles 
(open circles) and fast particles (closed circles) in the system as a function of the 
polymer volume fraction. 



to the incipient gel structure at surface sites but generally depart again before getting 
trapped. Preferential sites for this transient adsorbtion of fast particles are detected 
in pockets of the incipient gel formed by the slow ones, where they can establish more 
bonds with slow particles [37] . 

3.5. Implications for MCT 

It is interesting to point out that the average behaviour of the system is correctly 
described by MCT [33], both in the universal properties and the specific predictions 
for the attraction driven glass transition. The main exceptions to this involve the non- 
gaussian parameter (not shown; see [37]) and the discrepancy between 7 as measured 
from diffusion and from structural relaxation. However, the detailed dynamics, 
involving two concurrent populations of particles, signifies very strong dynamical 
heterogeneities which arc not really consistent with the approximation scheme lying 
behind MCT calculations. The strong DH observed in our incipient colloidal gels is 
related, it appears, to structural heterogeneities which have no direct analogue in the 
glass transition driven by repulsion [10]. At lower densities, these heterogeneities are 
even more important; in that region one expects to require significant modifications 
to MCT, as developed and discussed in Section 2. 
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4. MCT under shear 

The framework of MCT can be extended to address the relationship between stress 
and strain-rate in a system undergoing shear [9,38,39]; this includes a way to relate 
the yield stress of an arrested phase to its static structure factor S q as determined in an 
unstrained state. Via S q , the effects of attractive interactions on the nonlinear stress 
response can be incorporated. The response is always shear thinning. A schematic 
extension of this approach, allowing for a stress-dependent MCT vertex, admits the 
possibility of shear thickening states, which are often observed experimentally in 
dense repulsive colloids [40] . Although we expect attractive interactions also to have 
interesting effects on shear-thickening, we focus below on the shear-thinning case. 

The approach of [9, 38] predicts a finite yield stress a + at the arrest transition. 
For gels of attracting colloids, the scaling of a + with the width of the attraction 
is found below by combining a virial-type analysis of the AHSS model [41] with the 
model of 'isotropically sheared hard spheres' (ISHSM) [38, 39]. (This model isotropises 
advection in its effect on density fluctuations.) The resulting scaling can be compared 
to the value a + = 0.75fcsT/a 3 found at the glass transition of the ISHSM [39]. 

The ISHSM derives the nonlinear flow curves (viz., stress a versus strain rate 7) 
from the competition between the slowing down of the structural dynamics captured 
in the classical MCT, and the speeding up of fluctuations caused by shear advection. 
The latter mechanism has been intensively studied in the context of fluctuations in 
sheared systems close to criticality [42] , and close to ordering or microphase transitions 
[43]. It causes a time dependence of the wavevector of an arbitrary fluctuation 

q(i) = (q x ,q y + q x it,q z ) , (l) 

where the shear is along the x-direction with v x = jy; t here is the time since the start 
of shearing (or the birth of a fluctuation) . If a fluctuation initially has a wavelength 
~ 1/q, at a later time t its wavelength ~ |q(t)| _1 = l/q(t) will be smaller, leading to 
decay of the fluctuation caused by fast, small-scale particle rearrangements. 

The scaling of the yield stress a + with the attraction range S can, for small 

5, be estimated from a virial expansion of S q for the AHSS: one finds 1 — S^ 1 —> 
6A(j)S sm(2qa) / (qa) , where A = expe — 1. Inserting this expression into the ISHSM 
[39], one can take the limit 0^0 and A — > 00 with T v — 6(f>A 2 S/Tr 2 held fixed. This 
leads to the longitudinal memory kernel: 

mi(t) _ £ f' M (5_^ 2 cos p|zi) 4i(t)$K _ t| , (2) 

while the corresponding expression for the steady state shear stress becomes: 
. k B T 4>T V f°° 1 f k >~ -o (k{t)-k\~,^ 

«^i—y 1 * 1 ^ 2 cos(A^_j m . (3 ) 

Here, the limit S — > applies, and the rescaled wavectors q = qad and a cutoff fc> were 
introduced; choosing fc> w 3.68 maps the AHSS virial results onto the corresponding 
results of an attractive square well system with range <5 [41]. 

In the approach of Refs.[38, 39], these equations are closed via an MCT-like 
relationship for the normalized density correlation functions &q(t), whence we obtain 
predictions for the nonlinear rheology close to the arrest transition (at T° = 1.42) 
of attracting colloids with small 5 and low (f>. Although the low-0 approximation 
appears drastic, if shear is switched off in Eq.2, the results for the arrest line at small 
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5 agree qualitatively with those of Ref.[2] up to rather large 4> [12,41]. (Thus the 
approximations made are faithful to the standard MCT analysis, though they neglect 
the physics of heterogeneity discussed in Section 2.) With shear present, we can now 
analyse the altered behaviour of 3>q(i) at the onset of arrest and then compute the 
limiting stress at low shear rate by inserting the resulting asymptotics into Eq.3 [39]. 
This gives a finite yield stress a + just within the glass, which falls abruptly to zero 
on entering the fluid. 

Asymptotically the &q(t) obey a 'yielding scaling law', §q{t) — ► 3>J(i) with 
i = t/r{ A f) and T(j) a shear-rate dependent characteristic time. Also one finds the 
limiting closure relation [39]: 

= m+(t) | fdt' m+(t - $+(t) , (4) 

with m~~(i) obtained from Eq.2 using the limiting form (i). The yield stress er+ 
then follows from inserting <&~~(t) into Eq.3. 

The expressions for the memory kernel and stress, Eqs.2,3, contain a rapidly 
oscillating term cos(2(fe(t) — k)/8), which arises from interference of the particle 
density fluctuations within the narrow region of attraction. Studying this factor 
illuminates the role of shear. Without shear, constructive interference holds, and 
the factor is unity; the memory kernels describe as usual bond-formation owing to 
attraction, and the correlator stays arrested at its glass value, = $^(0) = /-. 

Under shear, the advection of wavevectors (giving k(t) — k = (k x k y /k)jt+. . .) produces 
rapid oscillations in this term when 6 is small. The interference is destroyed, causing 
a fast decay of the memory functions. The time scale r(7) needs to be found self- 
consistently from Eqs.2 and 4; the preceding argument shows that it scales with 
shear-rate and attraction range as r(j) — c a S /\j\ where c a is of order unity. In 
contrast, for repulsive interactions in the ISHSM, we found [38,39] r(j) — c r / 1 T I with 
c r again of order unity. Integrating Eq.3, one obtains (with further constants c' ar ) the 
scaling expressions c^iiss = c' a G a S and cr^g HSM = c' r G r for the yield stresses at the 
two transitions. Here G aj are the corresponding shear moduli which, within MCT, 
are predicted to acquire finite values on arrest. Their scalings have been discussed 
elsewhere [1,2,12,41,44]; for AHSS one has G a - kTa- 3 5- 2 whereas G r ~ feTa" 3 
holds for hard spheres. 

The final scaling result for AHSS is thus a + ~ G a S ~ fcTa -3 <5 -1 . This dependence 
of the yield stress on the attraction range leads to interesting scenarios in systems 
where both attraction and repulsion-driven arrest can be observed [7,45], and where 
at high enough concentrations the local structural dynamics should dominate the 
rheological behaviour. Interpreting the ratio of yield stress to the elastic modulus as a 
yield strain u y = er + /G, the scaling dependence of er + on 5 can easily be understood. 
The yield strain of a (high density) colloidal gel is of order the relative range of the 
attraction 5: the solid is shear-melted as soon as particle bonds are broken. At the 
same time, the scaling of the modulus is quite different for attractive and repulsive 
glasses; the attractive glass/gel is much stiffer as noted above. 

Thus two states of equal viscosity, close to the attractive and the repulsive 
branches of the arrest line, can have quite different relaxation times. It would be 
interesting to measure the dynamics, in the region where the two branches meet, 
along contours of equal zero-shear viscosity, as was done (without shear) for contours 
of equal diffusivity in recent simulations [46] . Because the yield strain is much smaller 
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for a bonded glass than for a caged one, the nonlinear rheology should vary strongly 
along these contours. 

A simple 'generalized non-linear Maxwell model' [39] summarizes the pertinent 
behaviour for both AHSS and ISHSM cases: 



Here 77(0) is the zero shear viscosity; 77(00) is a (small) limiting viscosity at very high 
shear rates (presumably set by hydrodynamic interactions, ignored here); and a* is 
a characteristic stress scale of order the yield stress discontinuity er + at the onset 
of arrest. The above model gives near-plastic yield in arrested systems (j](0) = 00) 
but can also describe a highly viscous but shear-thinning fluid phase [39]. In such 
a phase, close to a repulsive glass transition, the nonlinear viscosity 77(7) decreases 
significantly from its zero-shear value 77(0) once ijiO^a 3 /kT is of order unity. In a 
fluid close to an attraction-driven arrest, the decrease is only at much higher shear 
rates, with r}(Q) A /a? /kT w 1/8. This is associated with the much larger value of the 
yield stress (er + ~ 1/5), despite the much smaller value of the yield strain (u y ~ 8), 
in the nearby arrested phase. 

5. Concluding remarks 

We have described various developments in the application of MCT ideas to colloids 
with short range attractions. This theory is very successful at high densities in 
predicting re-entrant transitions induced by such attractions [7]; yet this re-entrance 
is caused by density fluctuations whose presence partly undermines the assumptions 
of homogeneity on which MCT is based. For a colloid volume fraction of <f> = 0.4 
we found by simulations that the MCT predictions for averaged properties remain 
generally good despite a surprising degree of dynamic heterogeneity. The latter is 
washed out in the averages, but can be probed more closely through the distribution of 
mean-squared displacements, and this analysis reveals a strong partitioning of particles 
into fast and slow populations. At much lower volume fractions the assumption of 
homogeneity becomes clearly wrong, but this can be partially addressed, at least 
at the level of qualitative prediction, by carefully combining MCT ideas with those 
of phase separation (arrest of one phase in a phase-separated morphology) and/or 
irreversible aggregation (applying MCT at the cluster scale). The effects of short 
range attractions on the rheology of colloidal suspensions are potentially quite subtle. 
The preliminary calculations reported above already show some interesting trends, 
particularly for the range-dependence of the yield stress in the attraction-driven glass, 
and for the shear-thinning behaviour in the nearby fluid phase. 
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